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DETERMINATION  OF  ELASTIC  MODULI  OF  SEA  ICE 

K.  L.  Williams,  R.  Stein,  T.  Wen,  R.  E.  Francois 

Applied  Physics  Laboratory 
College  of  Ocean  and  Fishery  Sciences 
University  of  Washington 
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ABSTRACT  in  conducting  experiments  tailored  to  the  needs  of  recently 

obtained  theoretical  results  for  longitudinal  and  shear  velocities 
For  applications  of  elastic  moduli  to  acoustic  prediction  models,  using  the  Biot  porous  media  theory  .20’'1  The  Biot  theory  is  a 

moduli  derived  from  acoustic  measurements  are  much  more  natural  one  for  porous  sea  ice.  Experimentally,  the  difficulty  lies  in 

appropriate  than  those  derived  from  mechanical  tests.  Several  the  need  to  specify  a  large  number  of  parameters  in  the  full  Biot 

researchers  have  measured  sound  velocity  in  ice  and  some  have  theory  for  an  anisotropic  open-pored  medium.22  However,  if  one 

used  these  measurements  to  deduce  ice  moduli.  Velocity  specializes  the  theory  to  examine  a  closed-pore  isotropic  porous 

measurement  techniques  include  resonance  vibration  of  ice  rods,  media  (which  is  probably  most  applicable  away  from  the  ice/water 

seismic  and  flexural  wave  measurement,  and  propagation  of  high  interface),  the  number  of  needed  parameters  reduces  to  seven,  all 

frequency  pulses  in  core  samples.  Differences  and  uncertainties  in  of  which  are  easily  interpreted  physically  and  are  amenable  to 

conclusions  to  be  drawn  from  these  studies  are  in  part  related  to  experimental  determination.  Even  with  this  specialization  the 

temperature  and  salinity  properties  of  the  ice  and  measurement  and  results  represent  a  generalization  of  previous  models.23  These 

core  handling  procedures  that  have  altered  the  structure  of  ice  from  theoretical  results  will  be  discussed  in  Section  III,  and  the  current 

its  in-situ  condition.  These  issues  indicate  the  advantages  of  a  status  of  the  knowledge  of  each  parameter  will  be  addressed, 

pulse  type  experiment  designed  to  allow  measurement  of  velocities  Conclusions  and  summarization  are  given  in  Section  IV. 
and  related  ice  properties  in  the  field.  We  pre'^nt  details  of  one 

such  design  including  the  equipment  used  and  the  resulting  moduli  II.  EXPERIMENT 
uncertainties  inherent  in  the  method.  The  theory  to  be  used  in 

interpreting  the  experimental  results  will  be  that  derived  by  Biot  The  primary  concern  in  the  experiment  is  measurement  of  acoustic 

for  a  porous  solid.  Expressions  for  the  longitudinal  and  shear  velocities  before  significant  temperature  change  occurs  in  the  ice 

velocities  in  ice  were  derived  assuming  a  sealed-pore  structure,  cores.  In  order  to  minimize  the  change  of  temperature  with  time, 

which  generally  applies  to  sea  ice  away  from  the  growing  the  cores  are  wrapped  in  closed-cell  foam  immediately  after 

ice/water  interface.  A  pilot  experiment  is  underway  at  the  present  withdrawal.  After  wrapping,  1  /8th  inch  diameter  holes  are  drilled 

time  to  perfect  the  procedure  and  will  allow  preliminary  into  the  core  at  regular  intervals  and  temperatures  recorded  with  a 

implementation  of  the  theory.  The  results  from  the  type  of  analysis  digital  thermometer  accurate  to  0.2°  C.  The  cores  are  then 

presented  here  ultimately  lead  to  prediction  of  the  acoustic  sectioned  in  a  mitre  box  (in  5-10  cm  lengths),  each  section  being 

impedance  layers  in  sea  ice,  a  crucial  factor  in  understanding  and  centered  around  the  temperature  monitoring  holes.  Starting  with 

predicting  acoustic  propagation  processes  involving  sea  ice.  the  section  furthest  from  ambient  temperature,  the  time  of  flight 

t  and  section  length  are  measured  using  the  apparatus  shown  in 

I.  INTRODUCTION  Fig.  1.  Identification  and  description  of  the  equipment  is  made  in 

the  figure  caption.  The  sending  transducer  at  one  end  of  the  sample 
Knowledge  of  the  vertical  profiles  of  acoustic  moduli  is  an  holder  is  excited  by  a  pulser/receiver,  and  the  delayed  signal 
important  prerequisite  to  such  tasks  as  ice  thickness  measurement  acquired  by  the  receiving  transducer  is  amplified  by  the  same 
or  determination  of  the  propagation  into  the  water  of  sound  pulser/receiver  and  output  on  a  10  MHz  digital  scope  where  it  is 

generated  on  the  ice  surface  .  It  is  also  implicitly  important  in  stored  on  the  built-in  disc  memory.  The  caliper  attached  to  the 

acoustic  studies  of  ice  keels  whose  raw  material  is  undeformed  sea  sample  holder  has  a  resolution  of  0.001  inches;  however, 

ice.  These  moduli  profiles  can  be  obtained  from  measurements  of  roughness  of  the  section  ends  and  use  of  a  coupling  gel  lead  us  to 

the  longitudinal  and  shear  acoustic  velocities  within  the  ice.  Three  use  a  more  conservative  estimate  of  length  uncertainty,  i.e., 
primary  velocity  measurement  techniques  have  been  used  in  both  measurements  to  0.01  inches  are  realizable  in  practice.  The  scope 
nonsaline  and  saline  ice:  seismic  and  flexural  wave  is  triggered  by  a  sync  pulse  sent  out  at  the  time  of  sender  excitation 

measurement,2-4  resonance  vibration  of  ice  rods,2,5-7  and  so  that  the  time  of  flight  is  directly  measurable.  Inherent  delay  and 

propagation  of  high  frequency  pulses  in  ice  cores. 7-1<>  "zero"  sample  length  are  checked  beforehand  with  sender  and 

receiver  in  direct  contact.  Also,  time  of  flight  measurements  on 
In  much  of  the  pulse  technique  work  ice  cores  were  taken  and  calibrated  aluminum  and  granite  rods  allow  elimination  of 
allowed  to  reach  equilibrium  with  the  air  temperature  or  were  systematic  time  scale  errors, 
refrigerated  and  taken  to  a  laboratory  for  the  velocity 

measurements.  The  multiphased  porous  nature  of  sea  ice  makes  it  After  the  time  of  flight  measurement,  the  temperature  of  the  core 

extremely  sensitive  to  this  kind  of  handling  procedure,  which  section  is  again  measured,  the  sample  density  is  determined,  and 

changes  the  in-situ  equilibrium  ratios  of  the  liquid  and  solid  the  sample  is  placed  in  a  sealable  container  and  allowed  to  melt, 

phases  as  well  as  the  properties  of  each  phase.  This  in  turn  changes  After  melting,  a  salinity  measurement  with  an  accuracy  better  than 

the  acoustic  behavior,  a  fact  that  provides  sufficient  motivation  for  0.1  ppt  is  made  using  a  conductivity  cell.24  Previous  APL  in-field 

design  of  an  acoustic  pulse  experiment  in  which  velocities  may  be  experimental  measurements  of  core  section  density  have  shown 

obtained  quickly  in  the  field  while  minimizing  and  monitoring  the  that  hydrostatic  measurements  of  densities  accurate  to  0.001 

change  of  temperature  profile  of  the  ice  core  and  minimizing  gm/cm3  are  achievable.  However,  the  present  method  is  time 

uncertainties  in  results.  In  Section  II  we  outline  the  experimental  consuming  and  there  is  some  small  unavoidable  temperature 

equipment  and  procedure  being  developed  and  tested  for  this  type  cycling.  The  technique  and  uncertainties  involved  are  discussed  in 

of  field  work.  We  also  address  briefly  the  uncertainties  in  elastic  Ref-  25.  The  main  concern  is  determination  of  whether  air  is 

moduli  that  result  from  experimental  limitations.  In  add-tm”  t~  entrapped  within  the  section.  Changes  in  density  due  to  air  content 

field  type  experiments  this  same  type  of  apparatus  will  be  useful  are  known  to  be  small  (on  the  order  of  0  03  gm/cm3);  however. 
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introduction  of  air  introduces  another  degree  of  freedom  that  may 
have  a  large,  frequency  dependent  effect  on  acoustic  velocities 
(see  the  end  of  Section  III). 

The  object  of  the  experiment  is  to  determine  ice  acoustic  moduli 
and  their  dependence  on  the  temperature  and  salinity  of  the  ice. 
Before  discussing  the  uncertainties  of  the  resulting  moduli  values, 
it  is  important  to  decide  which  moduli  are  most  appropriate  to 
determine.  Previous  analyses  have  concentrated  on  determining 
Young's  modulus  and  Poisson's  ratio.26  However,  for  comparison 
with  Biot  porous  solid  results,  it  is  more  natural  (see  Section  III)  to 
determine  the  effective  bulk  modulus  Keff1  and  shear  modulus  pc,<, 
of  the  ice  (w'hich  moduli  to  concentrate  on  is,  of  course,  simply  a 
matter  of  convenience  since  the  various  moduli  are  interrelated). 
These  moduli  may  be  determined  from  velocity  and  density  (p,) 
measurements  via  the  equations 


In  Section  III  these  effective  moduli  will  be  given  in  terms  of  the 
moduli  found  in  the  Biot  model  results  for  a  sealed-pore  porous 
solid.  Using  Eqs.  1  and  2  ,  the  measurement  uncertainties  given  in 
the  previous  paragraphs  as  the  limiting  errors,  and  assuming  an 
8  cm  core  section  length  with  q  =  4000  m/s  and  cs  =  2000  m/s,  the 
lowest  uncertainties  5pcffj ,  810^  obtainable  from  individual 
measurements  are  approximately 

8pcff,  =  1.3%  p'rf,  (3a) 

8Kcfr,  =  3.6%  Kcrf,  .  (3b) 

These  moduli  values  are  to  be  associated  with  temperature  (T)  and 
salinity  (S)  values  whose  uncertainties  are  also  given  in  the  above 
paragraphs  so  that  one  may  find  Keff,(T,S)  and  pcfrj(T,S).  Our  pilot 
experiments  are  aimed  at  perfecting  experimental  procedure  so  as 
to  approach  these  instrumentation  limits. 

Before  proceeding  to  theoretical  results,  it  is  appropriate  to 
examine  briefly  the  type  of  uncertainties  in  moduli  we  expect  due 
to  the  uncertainties  in  measurements  of  T  and  S.  This  calculation 
requires  knowledge  of  the  functional  dependence  of  the  moduli  on 
T  and  S,  which  is  actually  the  relation  being  sought  so  this 
uncertainty  determination  may  seem  somewhat  premature  (or 
impossible).  However,  we  can  obtain  an  estimate  for  the 
uncertainty  of  the  bulk  modulus  by  using  previous  results  for  the 
functional  dependence  of  the  Young's  modulus  of  sea  ice  on 
porosity26  ((3),  the  empirical  result  that  the  Poisson's  ratio  of  ice  is 
close  to  1/3  (Refs.  7,  26),  as  well  as  the  Frankenstein-Gamer 
equation  relating  P  to  T  and  S.  If  the  Poisson's  ratio  of  sea  ice  is 
assumed  equal  to  1/3,  then  the  effective  bulk  modulus  is  equal  to 
the  Young's  modulus.26  Using  Fig.  19  of  Ref.  26  to  determine  a 
linear  approximation  for  the  relation  between  porosity  and 
modulus,  then  relating  porosity  to  S  and  T  using  Eq.  13  of  Ref.  27, 
one  can  derive  the  relations 

Kcff,  =  l.OxlO11  -  (2.0xl0u)  p  (4a) 

K'ff,  =  l.OxlO11  -  (2.0xl09)  S  (  0.0532  -  .  (4b) 


In  this  relation  P  is  the  volume  of  brine  divided  by  total  volume  of 
an  ice  sample,  S  is  in  ppt,  the  Frankenstein-Gamer  equation  used  is 
appropriate  for  -0.5°C  <  T  <  -23°C,  and  the  resulting  modulus  is 
in  units  of  dynes  per  cm2.  The  approximation  in  Eq,  4a  is  limited 
to  values  of  p<0.5  since  the  bulk  modulus  must  be  positive.  The 
Bict  model  of  the  next  section  shows  that  the  P  range  of  this 
approximation  must  be  further  restricted  to  0<P<0.25  if  physical 
values  of  Biot  moduli  are  to  result  from  its  use.  Using  Eq  4b,  the 
uncertainty  in  Keffj  is  a  function  of  the  values  T  and  S. 
Uncertainties  for  a  few  combinations  of  T  and  S  are  shown  in 
Table  I. 


T(C) 

S  (ppt) 

P(%) 

6T  (C) 

6S  (ppt) 

6K/K  (%) 

-2.0 

4.0 

10 

0.2 

0.1 

3.1 

-2.0 

10.0 

25 

0.2 

0.1 

10  8 

-2.0 

10.0 

25 

0.2 

0.05 

1U.4 

-5.0 

4.0 

4 

0.2 

0.1 

0.5 

-5.0 

10.0 

10 

0.2 

0.1 

13 

-5.0 

10.0 

10 

0.2 

0.05 

1.1 

-20.0 

4.0 

1 

0.2 

0.1 

0.1 

-20.0 

10.0 

4 

0.2 

0.1 

0  1 

-20.0 

10.0 

4 

0.2 

005 

0.08 

Table  I  Uncertainties  in  the  effective  bulk  modulus  resulting  from 
uncertainty  in  measured  values  of  T  and  S  are  shown  for 
various  combinations  of  T,  S,  5S,  and  8T. 


If  the  time  of  flight  determination  of  the  bulk  modulus  were 
without  uncertainty,  an  individual  modulus  value  could  be 
associated  with  a  (T,S)  set  of  experimental  values  only  to  within 
the  uncertainties  given  in  the  last  column.  Thus  time  of  flight 
derived  modulus  uncertainties  (Eq.  3)  less  than  what  is  shown  in 
the  last  column  are  not  productive  until  equipment-limited 
uncertainties  in  T  and  S  are  lowered  The  uncertainty  in  Keff,  is 
largest  for  temperatures  near  freezing  and  high  salinities. 
Physically,  this  is  because  the  porosity  is  quite  large  (=  25  %)  and 
very  sensitive  to  temperature  changes.  Four  comments  are 
appropriate.  Eq.  4  from  which  the  table  is  derived  is  untested 
experimentally  at  porosities  above  about  10%;  experimental 
difficulties  due  to  brine  drainage  probably  represent  a  much  bigger 
challenge  at  temperatures  near  freezing  than  temperature 
uncertainty;  the  acoustic  velocity  determined  from  a  bulk  modulus 
known  to  within  10%  will  only  be  in  error  by  about  5%,  at  colder 
temperatures  time  of  flight  errors  are  the  major  contributors  to  bulk 
modulus  uncertainty. 

III.  BIOT  THEORY 

The  porous  nature  of  saline  ice27  is  responsible  for  its  differences 
in  behavior  relative  to  fresh  water  ice.  This  porous  nature  has  been 
taken  into  account  in  varying  degrees  in  some  analyses.  In 
particular  Refs.  15  and  23  discuss  calculations  of  sound  velocities 
in  which  the  saline  ice  is  modeled  as  a  matrix  of  solid  ice  with  the 
inherent  Young's  modulus  of  ice  and  pores  filled  with  brine.  A 
consequence  of  this  type  of  model  is  that  the  Young's  modulus 
depends  only  on  the  intrinsic  Young’s  modulus  of  ice,  the  bulk 
modulus  of  the  brine,  and  the  porosity  (Eq.  9.34  of  Ref.  23). 
However,  experimental  results  in  Ref.  7  for  Young’s  modulus  of 
summer  and  winter  ice  with  the  same  porosity  indicated  differing 
behavior  for  the  same  porosity.  This  was  tentatively  identified  as 
due  to  differences  in  the  solid  matrix  (thus  implying  use  of  the 
intrinsic  Young's  modulus  of  fresh  water  ice  for  the  matrix  might 
be  an  incomplete  representation). 

The  above  discussion  indicates  the  desirability  of  a  model  that 
allows  explicitly  for  the  variation  in  the  ice  matrix  as  conditions 
(i.e.,  T  and  S)  change.  The  Biot  model20'*1  of  a  porous  solid  is  one 
means  of  including  ice  matrix  variations.  However,  use  of  the  full 
Biot  theory  for  an  open-pored  anisotropic  material  requires  the 


2 


determination  of  17  parameters  22  At  least  initially,  it  is  useful  to 
specialize  the  most  general  Biot  results  to  reduce  the  number  of 
parameters  while  still  retaining  enough  complexity  that  the  results 
represent  a  generalization  of  previous  analyses. 


(8b) 


Reiterating  from  Section  I,  the  specialization  that  we  use  consists 
of  assuming  the  brine-filled  pores  in  the  saline  ice  are  sealed  and 
that  the  ice  is  isotropic.  The  sealed-pore  assumption  is  most  valid 
away  from  the  skeletal  layer  at  the  water/ice  interface.27  With  this 
restriction  the  number  of  parameters  that  must  be  specified  is 
reduced  to  seven,  which  are  easily  interpreted.  Before  synopsizing 
the  results  derived  using  Biot  theory,  it  is  appropriate  to  point  out 
that  since  the  model,  as  derived  by  Biot,  was  for  a  tw'o- phase 
problem  (i.e.,  liquid/solid,  gas/solid,  etc.)  it  may  not  reliably 
predict  behavior  in  a  three-phase  situation  (liquid/gas/solid). 
Further  comment  about  the  three  phase  situation  will  be  made  at 
the  end  of  this  section. 

The  terminology  to  be  used  here  mirrors  that  of  Ref.  28.  A 
complete  development  of  Biot’s  theory  may  be  found  in  Refs  20, 
21,  29,  and  30.  We  simply  note  here  that  Biot's  examination  of 
propagation  of  energy  through  a  porous  material  began  with  a 
determination  of  the  potential  energy  stored  within  the  material 
and  the  kinetic  energy  due  to  the  motion  of  both  the  solid  and 
liquid  phases.  With  these  he  could  write  the  Lagrangian  density 
and  determine  the  equations  of  motion.  From  this  type  of  analysis 
one  can  derive  a  set  of  coupled  equations  of  motion  for  the  liquid 
and  solid.  These  equations  are  given  as  Eq.  6  of  Ref.  28  (we  note  a 
typographical  error  in  the  first  of  those  equations,  i.e.,  2p  should  be 
replaced  by  p).  These  equations  are  coupled  via  the  relative 
motion  of  the  liquid  and  solid.  Assuming  a  sealed-pore  structure 
eliminates  this  relative  motion.  To  correctly  handle  this  situation 
within  the  context  of  Biot  theory,  one  must  go  back  to  the  potential 
and  kinetic  energy  equations29  Wand  eliminate  the  relative  motion. 
In  doing  so  one  finds  that  only  one  equation  survives  and  is  of  the 
form 


pb  V2!?  +  (H-pb)V (Vif)  =  p  ,  (5) 

8‘t 


where 


P  =  PPf  +  d-p)Ps 


Substitution  of  the  expression  for  H  (Ref.  30)  gives  our  desired 
results 


/  |(Kr-Kb)2/(Kr(l  +  p(K/K(  -  1))-Kb)|  Kb  +  4pb/3 
V  ppf+(l -P)ps 


In  this  expression,  in  addition  to  previously  defined  quantities,  Kr 
is  the  bulk  modulus  of  nonsaline  ice,  Kf  is  the  bulk  modulus  of  the 
pore  brine,  and  Kb  is  the  bulk  modulus  of  the  ice  frame.  These 
expressions  are  novel  in  that  they  allow  one  to  explicitly  see  the 
dependence  of  the  acoustic  velocities  on  porosity  and  frame 
rigidity.  Comparison  of  Eqs.  9  and  10  with  Eqs  1  and  2  snows  thai 

Pcff,  =  Pb  01) 


'Cff  _ 


(Kr-Kb)2 


(Kr(l  +  P0VK,-  1))-Kb) 


+  Kb 


p,  =  Ppf+(l  -P)ps  .  (13) 


The  results  in  Eqs.  9-13  represent  a  generalization  of  results  based 
on  simpler  theories.  In  particular,  to  recapture  velocity  results  from 
suspension  theory32  one  sets  Kb  and  pb  equal  to  zero.  To  obtain  the 
Young’s  modulus  in  Eq.  9.34  of  Ref.  23  (see  beginning  of  this 
section),  one  assumes  Poisson’s  ratio  is  zero23  and  sets 
Kb  =  (l~P)Kr  .  Setting  Kb=  (l-p)Kr  is  equivalent  to  assuming  the 
ice  wi'hfn  th'  frame  matrix  has  a  bulk  modulus  equal  to  the 
inhere’  .n^jlus  of  ice.  Clay  and  Medwin33  have  derived  an 
express;,  cv  alent  to  Eq.  12  via  a  different  method. 


In  these  equations  u  is  a  displacement  vector  for  differential 
volume  of  the  solid,  pb  is  the  shear  modulus  of  the  frame,  p  is  the 
total  mass  density,  pf  is  the  fluid  density.  p5  is  the  solid  density, 
and  H  is  the  combination  of  bulk  moduli  for  the  frame,  solid,  and 
fluid  given  below.  If  the  displacement  vector  is  written  as  a 
combination  of  scalar  (<»  and  vector  (y)  potentials 


u  =  V<$  +  Vxy  , 

(6) 

then  Eq.  6  can  be  rewritten  as  two  equations 

HV20  =  pO 

(7a) 

pbV>  =  py 

(7b) 

where  the  double  dots  represent  partial  differentiation  with  respect 
to  time.  Assuming  monochromatic  solutions  for  the  potentials31 
one  can  arrive  at  the  following  longitudinal  and  shear  wave 
velocities 


The  results  aoove  indicate  that  seven  parameters  ate  needed  to 
characterize  a  closed-pore,  air-free  Biot  solid,  i.e., 
P,  ps,  pf,  pb,  Kb,  Ks,  Kf.  All  are  easily  interpreted  physically.  Of 
the  seven  parameters,  two  depend  on  temperature  (Ks,ps)  and  the 
rest  on  both  salinity  and  temperature.  It  is  anticipated  that 
characterization  of  the  ice  cores  examined  via  the  experimental 
procedure  of  Section  II  will  be  completed  as  follows.  The  values 
of  S  and  T  will  be  used  to  determine  p,  pf  and  ps(T  only)  using 
Ref.  34.  The  consistency  of  these  predictions  with  the 
experimentally  determined  ice  density  will  be  checked  using 
Eq.  13.  Low  experimental  values  for  ice  density  imply  air  content, 
and  sections  containing  air  (most  likely  in  the  upper  few 
centimeters  of  the  core)  will  not  be  analyzed  in  terms  of  the 
present  theory  (further  comment  on  the  modeling  of  a 
liquid/air/solid  porous  medium  is  made  at  the  end  of  this  section). 
For  air-free  sections  the  experimentally  determined  effective 
moduli  will  be  used  to  determine  pb,  and  Kb.  Obtaining  pb  is 
trivial  (cf.  Eq.  11),  but  determination  of  Kb  requires  some 
explanation.  First,  Eq.  12  can  be  solved  for  Kb  to  give 


Kb  =  Ks 


pKctt,(Ks-Kf)  -  K,<KS-Kcfr,) 
pKs(Ks-Kf)  -  KffKs-K'",) 


(14) 
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Equation  14  shows  that  in  addition  to  the  experimentally 
determined  value  of  Keffi(S,T),  we  need  KS(T),  and  K,<S,T)  to 
obtain  Kb(S,T).  Measurements  of  the  temperature  dependence  of 
Ks  may  be  found  in  Refs.  5  and  8.  No  measurements  of  the  salinity 
and  temperature  dependence  of  Kf  for  salt  water  brine  have  been 
found,  and  at  present  we  will  assume  Kf=2.0xl010  dynes/cm2. 
This  value  corresponds  to  an  acoustic  velocity  in  water  of  about 
1440  m/s.  Measurement  of  both  KS(T)  and  Kf(S,T)  are  needed  and 
are  amenable  to  well  controlled  laboratory  experiments,  implying 
experimental  errors  much  lower  than  those  inherent  in  the  current 
measurement  of  effective  moduli.  A  preliminary  estimate  of  the 
error  in  the  values  of  Kb  due  to  the  errors  in  Kcff,  can  be 
established  using  Eq.  14,  the  approximation  of  Eq.  4a,  and  using 
the  fact  that  Kf  is  found  to  equal  approximately  0.2  Ks.  With  these 
assumptions  Kb  =  (1-4(3)  Ks  and  the  partial  derivative  of  Kb  with 
respect  to  K  j  is  found  to  equal  approximately  4,  implying  an 
error  in  measurement  of  Kb(S,T)  of  0.14  Kcffj(S,T)  for  the  error  in 
Eq.  3.  The  result  for  Kb  implies  that  the  relation  in  Eq.  4a  cannot 
be  used  for  (3>0.25  since  physically  0  <  Kb  <KS.  It  is  important  to 
emphasize  that  one  of  the  major  uses  of  the  Biot  moduli  will  be  to 
predict  acoustic  velocities  in  ice,  and  if  the  error  in  Kb  is  traced 
backward,  it  implies  velocity  results  to  less  than  2%  error. 

A  comment  on  inclusion  of  air  in  the  porous  solid  model  is  in 
order  to  conclude  this  section.  Bedford  and  Stem35  have  extended 
the  Biot  model  to  include  content  of  a  small  volume  of  air  bubbles. 
Their  results  indicate  that  above  the  resonance  frequencies  of  the 
air  bubbles,  the  velocity  of  the  longitudinal  wave  examined  here 
would  be  essentially  the  same.  However,  below  the  bubble 
resonance  frequency,  the  velocity  would  be  substantially  depressed 
and  at  resonance  the  velocity  goes  through  a  maximum.  The 
applicability  of  their  analysis  in  obtaining  ice  acoustic  velocities  is 
currently  being  examined. 

IV.  SUMMARY  AND  CONCLUSIONS 

An  experiment  currently  being  tested  for  in-field  measurement  of 
ice  moduli  has  been  outlined.  The  in-field  nature  of  the  experiment 
is  a  requirement  of  the  sensitivity  of  saline  ice  to  handling 
procedures.  The  moduli  to  be  extracted  are  those  needed  in  the 
Biot  model  specialized  to  an  isotropic,  sealed-pore,  air-free  porous 
solid.  The  specialization  reduces  the  number  of  parameters  needed 
to  seven  while  still  generalizing  previous  results.  Operationally, 
the  end  product  of  these  kind  of  experiments  is  an  increase  in  the 
accuracy  of  measuring  and  predicting  velocity  profiles  from 
relatively  simple  measurements  of  the  physical  parameters  (S.T), 
and  therefore  an  increase  in  the  accuracy  of  analyses  and  tasks  that 
use  them. 

The  Biot  model  has  been  experiencing  increased  use  in  recent 
years.  Its  major  drawback  is  the  need  to  experimentally  measure 
many  parameters.  Since  the  model  is  only  as  good  as  the  parameter 
measurements,  it  is  important  that  experimental  parameter 
measurement  efforts  be  initiated.  The  work  here  is  in  part 
motivated  by  the  need  for  accurate  parameter  determination. 
Experience  with  the  present  model,  in  which  the  number  of 
parameters  are  restricted  and  easy  to  interpret  physically,  is  an 
important  first  step  to  measurements  of  parameters  in  more  general 
situations. 
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Figure  1.  The  in-field  ice  velocity  measuring  apparatus:  a)  left  lo  right  -  Nicolet  4094  digitizing  scope  with  a  10  MHz  plug-in 
unit  and  built-in  disc  memory.  Panametrics  5055PR  pulser/rccciver  (p/r),  APL  built  ice  sample  holder  that  houses  sending  and 
receiving  transduccrs;b)  typical  received  waveform  after  amplilication  by  p/r  unit;c)  close-up  of  sample  holder.  In  operation  the 
p/r  excites  a  broad  bandwidth  transducer  (longitudinal  or  shear)  coupled  to  the  ice  sample.  The  receiver  signal  is  amplified  by  the 
p/r  and  the  resulting  signal  displayed  on  the  scope.  Sample  length  is  directly  measurable  by  the  caliper  attached  to  the  sample 
holder  and  time  of  flight  is  dirccdy  measurable  on  the  scope.  The  emphasis  is  on  fast/accurate  data  acquisition  so  that  the  physical 
parameters  of  the  sample  are  near  their  in-situ  values. 
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